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In this paper, the product of parabolic cylinder functions D v (fiiz)D v + tl -\{z), with different 
parameters g and u, are established in terms of Laplace and Fourier transforms of Rummer’s 
confluent hypergeometric functions. The provided integral representations are transformed 
to easily yield Nicholson-type integral forms and used to derive other series expansions for 
products of parabolic cylinder functions. 

Keywords: Parabolic cylinder functions; Weber-Hermite’s functions; Confluent 
hypergeometric functions; Laplace and Fourier transform; Nicholson-type integrals. 

AMS Subject Classification: 33C15; 44A10; 44A20; 30D10; 30E20. 


1. Introduction 


Parabolic cylinder function, called also Weber-Hermite’s function [1], arises in different 
scientific disciplines and is considered in many branches of mathematics. For instance, 
in number theory, Helfgott has recently utilized it as a tool in his proof of Goldbach’s 
weak conjecture [2], He developed a new asymptotic expansion of the parabolic cylinder 
function using saddle point method [3, 4]. It is also widely considered in physics thanks 
to its strong relation with the Gaussian normal distribution. Actually, it represents the 
general moments of the normal distribution; see its integral representations in [5, p. 1028]. 
Parabolic cylinder function is the solution of the differential equation 


fy, 

dz 2 



(1) 


and is entire function of all z 6 C and v G C. 

Following the notation in [5-10], we use D u (z) to denote the parabolic cylinder function. 
Note that other notations exist in literature such as U(u,z ) = D_ u _ l / 2 (z) [2-4, 11— 
13]. It is known that, for nonpositive integer v , D v (z) and D u (—z) are two independent 
solutions of the differential equation (1). Also, D u (z) and D- v -i(iz) are two independent 
solutions for any v € C [6, p. 117, Eq. (10), (11)]. 

Parabolic cylinder function D u (z) may be expressed in terms of Hummer’s confluent 
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hypergeometric functions <h and T by 
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A number of integral representations and series expansions of parabolic cylinder func¬ 
tions are given in [4-6, 11-13]. Moreover, integral representations of the product of two 
parabolic cylinder functions have been also investigated and several formulas are provided 
in [5, 6, 14]. Among them, it is known that (see [5, p. 697] or [6, p. 120]) 



where J u is the Bessel’s function of the first kind. Other integral representations of the 
products of D u , involving modified Bessel’s function of the second kind, are also found 
in [5, 6]. 

Furthermore, Durand in [15] derived a Nicholson-type integral for the product of two 
parabolic cylinder functions, valid for > — 1 and for all zG C. 
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Malyshev provided in [7] other useful integral representations for the products 
D u (z)D u (±z), valid for < 0. He showed that 


D v (z)D v (±z) 
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It is important to notice that the integral representation (5) of D u (z)D u (—z) provided 
by Malyshev is equivalent to Durand formula (4) because, using the relation between 
D u (z) and D_\- V (iz) in equation (24), it is easy to show that 


D u (z) 2 + —J — (cos ttuD u (z) - D u (-z)) 2 = + D-i_ u (iz)D~i_ u (-iz) (6) 
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Therefore, replacing u by — 1 — v and z by iz in (4) leads to (5) for D u (z)D I/ (—z). 


Very recently, Glasser provided in [8] an integral representation for the parabolic cylin¬ 
der function product D u (x)D u (—y) with unrelated arguments x and y satisfying x > y. 
In the same direction, Veestraeten [9] has recently derived other integral representations 
for D v (x)D u {y) and D v {x)D v -\(y) with unrelated real arguments x and y satisfying 
x + y > 0 . He has also pointed out that, using the recurrence relation of the parabolic 
cylinder function, it is possible to extend his results to have the same kind of integral 
representations for D u (x)D u _ n (y), where n is an integer. 
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Furthermore, Nicholson-type integral (4) is useful in the study of real zeros of parabolic 
cylinder functions. In their work [1] about the properties of real zeros of any linear com¬ 
bination of solutions of (1), Elbert and Muldoon established a remarkable identity that 
can be written, after a straightforward transformation, in terms of parabolic cylinder 
functions as follows 
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where erf(z) = 1 — ^2/ite~ z2 !' 2 D-i [z\f 2) is the error function. 

Identity (7), valid for Sfiz/ > — 1, was the key tool to obtain in [1] an explicit formula for the 
derivative of a zero of Hermite’s function with respect to the parameter. Besides, as can 
be observed, identities (4), (5) and (7) seem to be special cases of a general formula. To 
the best of my knowledge, no such general identity for the product of parabolic cylinder 
functions with different parameters exists in literature. 

Primarily motivated by the need to generalize identities (3) to (7), we have in this paper 
three aims. The first is to establish general formulas covering in particular identity (3). 
As will be shown later, the new formulas explicitly give the product of parabolic cylinder 
functions with different parameters and involve the Laplace and Fourier transforms of 
confluent hypergeometric functions. The second aim of this work is to transform the new 
found identities to yield general Nicholson-type integrals for the product of parabolic 
cylinder functions with different parameters. These identities are the generalizations of 
(4), (5) and (7). In addition, we extend the results from the first and second aims to 
point out some series expansions of the product of parabolic cylinder functions. 


2. Integral representations of products of parabolic cylinder functions 

In what follows, we establish the Laplace and Fourier transforms of confluent hyperge¬ 
ometric functions of argument x 2 . We show that these integral transforms are explicit 
and equal to the product of parabolic cylinder functions with different parameters. More¬ 
over, general Nicholson-type integrals for the parabolic cylinder functions products are 
deduced. 

Throughout the paper, we denote by and T the Rummer’s confluent hypergeometric 
functions of the first and second kind, respectively. 

Let [y) n = T{y+n)/T{y) be the Pochhammer symbol, confluent hypergeometric function 
<F is analytic on C and is defined, as in [5], by the series expansion 
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Equation (8) is equivalent to the following integral representation of $ 


$(i/, n- z ) = — -- /’Vv'- 1 (l - tY~ v ~ l dt , 0 < Tiy < T/x 

r w)r (fi-v) Jo 


( 9 ) 


3 



On the other hand, confluent hyper geometric function of the second kind T is defined as 
in [5, 6] by 

M z) = ^ 1 ^ $(^, n\ z ) + r ^ I/ ) 1 ^ 1 ~ A * $ ( I/ + 1 - /b 2 - //;*) (10) 

1 /-+0O 

= —- / , ?ftz/ > 0, ftz > 0 (11) 
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2.1. Laplace and Fourier transforms of confluent hypergeometric functions 

In this subsection, we derive some useful identities of the product of parabolic cylinder 
functions involving Laplace and Fourier transforms of confluent hypergeometric func¬ 
tions. 


Theorem 1 Let D v and be defined as in (2) and (8), respectively. Let a 0 be a 
purely imaginary complex number, i.e., 9ia = 0. Let v and p be two complex numbers 
such that iRp > 0. Then, for every complex number z satisfying iRz > 0, we have 

1 r+°° i , 
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Proof. Let v and p be two complex numbers such that 0 < IRu < tkp. It follows from the 
integral representation (9) of that we have, for every nonnegative real number x and 
purely imaginary complex number a, 
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The passage to the last integral is just a matter of change of variable, i.e., xt replaced 
by t. 

By writing x 2 — 2 xt = (x — t) 2 — t 2 , it is immediate that 
—^a^e-S 0 * 2 * ( v,p;ax 2 

=f * g{x) (14) 


) =— - y -- l\\ at ^- t ^- 1 (x-tY-^dt 

T(u)T(p-u) J 0 


where f * g is the Laplace convolution of functions / and g, both defined on nonnegative 
real number t by respectively: 

aP at 2 

m = wf 2 and 3{t) = nPPf 2 ' 
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It is known that the Laplace transform of a convolution integral is the product of the 
Laplace transforms of the two functions [5, p. 1108]. 

Let &(z) be the integral in the right hand side of (12). It then follows from (14) that 

Hz) = £f(z)Cg{z) (15) 

where Cf(z) is the Laplace transform of function /. 

From the integral representation of parabolic cylinder functions in [6, p. 119, Eq.(3)], we 
have 


£/(*) 
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and 


£g(z) = a 2 2 e L D v -n 
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(17) 


Making the product of (16) and (17), we complete the proof. ■ 

Note that the restriction made on v and p at the beginning of the proof is required 
only for the integral representation of confluent hypergeometric function 4>. There is no 
need for any restriction on u in identity (12), it is valid for any v 6 C and p € C such 
that IRp > 0. The restriction made on a is required because the integral in the right hand 
side of (12) diverges if SRa 0. 

Theorem 2 Let D v and T be defined as in (2) and (10), respectively. Let a 0 be a 
complex number with Ka > 0. Let v and ji be two complex numbers such that 0 < < 2. 

Define '(/(z) be the following integral transform 
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Then for every complex number z, we have 
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and consequently 


D- v 
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Proof. To prove (19), we assume first that $ta = 0 and SRz > 0. 

Starting with the representation (10) of 'll in terms of and appealing the relation 
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r(/i)r(l — /j) = ttcsctt/j,, we have the following computation. 
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Using the result of theorem 1, it directly follows that 
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The integral (18) of 4/(2:) is convergent for every complex a / 0, such that ffia > 0. 
Then, by analytic continuation, identity (22) holds also for any complex a /0 satisfying 
> 0 because the right hand side of (22) is valid for every complex number a / 0. 
Moreover, it is worth noting that when = 0, the restriction on ^Rz > 0 is required 
for the convergence of the integral in the right hand side of (18). It is also clear that the 
integral in (18) is convergent only for 0 < Sfi/x < 2. 

Now, multiplying both sides by -— s ‘ n7rM and using the constraint Ka > 0 yield 


a 2 sin 7174 
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For real parameters an arguments a and 2 :, the left hand side of (23) is real so also the 
right hand side. Then, replacing the right hand side of (23) by its complex conjugate is 
also valid. 

The proof of identity (19) results directly from (23). We have to just take the right hand 
side of (23) and make some transformations using the following known relation between 
D_ u and A-i [6, p. 135]. 


e" ( " 1) D_ u (iz) + e l) D- v (-iz) = f^j D „-i(z) 


(24) 


Then, it follows from (23) that 
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Simplifying (25), we arrive at identity (19). 


To prove (20) from (19), we recall the integral representation (18) of T(z), substitute 
z with —iz, equate from each side even parts on z together and then odd parts also 
together, we find after simplification the following Fourier Cosine and Sine transforms: 

Li +oo 

sin f cos(zx)x fJ, ~ 1 e~ 2 ax T (y, n\ ax 2 ) dx 

= D - (*) D ^- (-*) + D - (-*)(*) < 26 > 

and 

An 2 /*+oo ^ 

^— cos ^ j sm(zx)x fJ- ~ 1 e~ 2 ax T (i/, /a] ax 2 ) dx 

= D - (*) (-*) - *>-, (-*) *v>- (*) ( 27 > 

Adding (26) to (27) and using that sin ^ cos(za:) + cos ^ sin(zx) = singer + ^), we 
arrive at (20). 

It is important to note that in contrast to (18) and (26) where [i should satisfy 0 < Sfi/i < 
2, the left hand side of (27) converges for — 1 < Sfi/i <3. ■ 

Recall form [5, p. 1028] that Bessel’s functions are special cases of confluent hyperge¬ 
ometric functions. 




Jv- 1/2(t) = " ~ le ^ ® (v^ix 2 ) 

1 [V + 2 J 


(28) 


K U- 1 / 2 ^) = Vnx 2u l e {v, 2V]X 2 ) 

Set then n = 2u, a = i in (12) and a = 1 in (20) and use T duplication formula 


o2z/-l 

T{2 V ) = — f =-T{y)T{y+\) 


(29) 


(30) 


we have immediately the following identities, already known in literature [5, 6]. 

D_ v (^zeD-„ e ~ ZXj u _l (f) dx , Ki/ > 0, ^z > 0 (31) 

and 

/o r+oo , .. 

D- u (z)D u _i(—z) =— / sin (zx + ttv) K i ( ) dx , 0 < < 1 (32) 
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Remark 1 The integral transform 'P(z) in (18) can be analytically continued for all 
R/a > 0 by just removing the singularity of T(zv, /i; x 2 ) at x = 0 as follows (for a = 1). 


T(z) 



^x M 1 'h [y, y x 2 ) 


- !) (y + 1 - h)k x 2k+1 A 
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where n is any natural integer that should satisfy n > L^r^J • 

Theorem 3 Let D u , <J> and T be defined as in (2), (8) and (10), respectively. Let v be a 
complex number such that Ru > 0. Then, the following integral transforms hold for every 
complex number z. 

z 2 o 2v r+oo 

eTD^(z) 2 = -4— / ^-i e -2—2* 2 * (l/>1/ + i ;x 2) dx (34) 
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Proof. To prove (34), let Kz/ > 0 and observe that the function x i->- 

x 2u ~ 1 e~ 2x ”4? (y, v + ),; x 2 ) can be seen as a convolution integral. We have actually form 
the integral representation (9) of <J> that 
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Making the change of variable u : t >-+ x\/l — t yields 
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Applying now Laplace transform to the function x [f * /]( 2x) and appealing the 
integral representation of parabolic cylinder function from [6, p. 119, Eq.(3)], we obtain 

vrrvw / e~ 2zx [f * f]{2x)dx = e^D_ v (z ) 2 (39) 

1 W Jo 

Using the duplication formula (30) of T function, we complete the proof of identity (34). 

To prove (35), we use the following identity: 

Let g and h be two summable functions on [0, +oo), then we have 

/*+oo /•+oo z'+oo /»+oo /*+oo /» +oo 

/ g{x)dx / h(x)dx = / g(x)h(x + t)dxdt + / / h(x)g(x + t)dxdt 

Jo Jo Jo Jo Jo Jo 

(4 °) 

Its proof is simple. Consider the primitive integral functions G and H of respectively g 
and h: 


G{x) 


r*+oo 


g(t)dt , and H(x) 



h(t)dt 


It is clear that G(+oo) = H(+ oo) = 0. Also, observe that the right hand side of (40) can 
be seen as 


r-+oo 


G\x)H( x) + G(x)H'(x)^ dx = G(0)H(0) 


Applying now relation (40) to g(t) = t v l e zt 2 and h(t) = t v l e zt 2 yields 
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Changing the variable x with — (\Ju + 1 — 1) in the inside integral of (41), we have 

dx = -rr =r and thus 
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Changing the variable t with 2f in (41) and using (42), we complete the proof of (35). 

Identity (36) is a direct consequence of (34) and (35). To prove it, we shall use the 
relation between T and $ in (10) and apply some transformations utilizing identity (24). 
We do not need any restriction on v for identity (36); it is valid for every v £ C. ■ 

I tried to generalize (34) to (36) for the product of parabolic cylinder functions with 
different parameters v and g but it seems to do not involve integral transform of confluent 
hypergeometric functions. 


9 



2.2. Nicholson-type integrals for products of parabolic cylinder functions 

As explained in the introduction of this paper, we interest to generalize identities (4) 
to (7). The idea is to transform the results from the previous theorems to obtain the 
following corollaries. 

Corollary 1 Let D u be defined as in (2). Let // and u be two complex numbers such 
that Ru < 0. Then we have for \arg(z)\ < \ 


D v {z)D v ^ x {z) = -^1 r°° e (^+u)t-^ cotht^_ i f zVco ^\ 

0 v 7 


dt. 


(sinh 2 1) 2 


(43) 


Proof. Here, we employ again the result of theorem 1 and we recall the integral represen¬ 
tation of confluent hypergeometric function in (9). Assume temporarily that 0 < v < g, 
a € iR* and z > 0, we have thus the following direct computation. 


v u 
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(44) 


The last equation is obtained by firstly interchanging the order of integration in the 
double integral because both converge for the condition made on v, g, a and z\ and sec¬ 
ondly by evaluating the inside one using the integral representation of parabolic cylinder 
function from [6, p. 119, Eq.(3)]. 

We observe also that the last integral in the right hand side of (44) is convergent for a 
not necessary purely imaginary complex number assuming of course some restrictions on 
z. Hence, by analytic continuation, the right hand side of (44) is also equal to the left 
hand side for positive real number a. This enables us to set a = 1 in what it follows for 
the proof. 

Making the change of variable x = 1 — 2t and splitting the integral into two parts, we 
have 


Tciv TV n\ 2 1 u r i 

D_ u (iz)D u ^(z) = -r / e^H_, 
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Making once again the change of variable x = tanh t in both integrals of the right hand 
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side of (45) and simplifying, we arrive at 
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We take imaginary parts from each side of (46) using identity (24). It follows that 

2 1 2 p+°° + „ / t —;—t\ dt 


D u _i(z)D u _ ll (z) = -v) J e {2l/ ^ 4 cotht D^i (^zYcotht^j 


(sinh 2t) 2 


(47) 


Replacing u by v + y gives the desired identity (43). 

The integral of (43) is convergent at t = +oo for Rxx < 0, while at t = 0, it converges 
provided that \arg(z)\ < \ or R/x < 2 and \arg(z)\ = f • ■ 


Equation (43) is the generalization of identity (5) because with /x = 1 in (43), we read¬ 
ily find the Nicholson-type integral of D u (z) 2 which was derived by Malyshev in [7, p. 7, 
Eq. (23)] and extended by Glasser in [8] or by Veestraeten in [9] for unrelated arguments. 

In addition, we extend results of theorem 2 to establish the following corollary that 
generalizes identities (4) and (7). 

Corollary 2 Let D u be defined as in (2). Let /x and u be two complex numbers such 
that R/x < 2 and Rzx < 0. Then, for every complex number z, we have 

= e^X-WTV. (*v/taSi) (48) 


Proof. We proceed as in (43). To prove (48), we use the result of theorem 2 with a = 1 
and we recall the integral representation of confluent hypergeometric function T in (11). 
It follows then from (20) that 


D_ u (z)Dn_ u _i(-z) 

2 /»+oo /*+oo 


\/ 2 nT{y) Jo Jo 


n-\-oo /»-hoo ^ 

/ / sin (zx + x»- 1 e-2 x2 ( 1+2t k l '- 1 (l + tY~ v ~ l 

Jo Jo 


dtdx 


(49) 


Interchanging the order of integration in (49), evaluating the inside integral (the one 
with variable x) using the integral representation of Z) /i _i(— z) in [6, p. 120, Eq.(4)], 
then making the change of variable y = and simplifying, we obtain 


D_ u {z)D^- v -\{-z) 


2^-n r 1 2 

W)L e ' D "- i( - z ^ ) 
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(50) 
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Changing the variable of the integration y = tanhf as in (45) and simplifying, we have 


D-vWDp-v-^-z) = ^4 tanht^_ i ^ (51) 

1 (Rj Jo v 7 (smh2t) 2 

Replacing u by — u and z by —z, we complete the proof. 

The integral of (48) is convergent at t = +oo for $lu < 0, while at t = 0, it converges 
provided that $ty <2. ■ 

Remark 2 Using the differential equation (1) satisfied by the parabolic cylinder function 
and making the first four derivatives of the function y : z i—> D u (z)D u+fl -i(z) with 
respect to z, it is easy to prove that the function y satisfies the following fourth order 
differential equation 

g + 4(„ + S ^)g-34 +M ,-2)^0 (52) 

An alternative proof of corollary 1 is then to show that the right hand side of (43) 
respects also the differential equation (52) and both sides of (43) are equal for the initial 
conditions, i.e., y^(0), where 0 < n < 3. 

Similarly, one can prove that both sides of equation (48) satisfy also the fourth order 
differential equation (52) and are equal at the initial conditions. 

Now, setting y = 1 in (48) directly provides Nicholson-type integral (5) of D u (z)D u (—z) 
that was established by Malyshev in [7, p. 7, Eq. (26)] and by Durand in [15] with another 
form as explained in the introduction. 

It would be nice to outline other particular cases of identity (48). Then, if we take from 
(48) the even and odd parts separately using the expression of D u in (2), we establish 
the following Nicholson-type integrals: 


D u (z)D l/+ll _i(-z) + D u (-z)D u+fl _ 1 (z) 
2v / 2t r 


f+°° e (2iH-/d* 


r(-i/)T(l-f)y 0 (sinh2f)§ 


4> 


( 2,2; ~ 1> 2 tanli 4 


dt 


(53) 


D„(z)D u+tJ ,-i(-z) - D u (-z)D u+fJi _i(z) 
4 Zy/^ [+°° e ( 2 ^+^)i 


T(-z/)T(±- %)J 0 (sinh2f)§ 


V tanh + T t an h 1) 


dt 


(54) 


Equation (53) is valid for < 2 and $tu < 0, while (54) is valid provided that Ry < 3 
and 'Sbs < 0. However, it is possible to remove the restriction on y for (53) and (54). 

We take the derivative of both sides of (54) with respect to y and then set y = 1, we find 


D„{z) 


dD v {—z) 

du 


D u (-z) 


dD v {z) 

du 


( 1 3 . 

T(—u)J 0 coshf V ’ 2 ’ 


z ? 2 tanh t\ dt (55) 


where -^=<3? (l, — z 2 ) can be reduced to ie z2 erf(iz); with erf(z) is recalled the error 

function. 

Identity (55), valid for ffiu < 0, is equivalent to [1, p. 8, Eq.(5.13)] provided by Elbert 
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and Muldoon. 

Similarly, the derivative of (53) with respect to fi immediately leads, for // = 2 after 
removing singularity, to 




y/2i tz 2 


r-+oo g 2(z/+l)i 


dv 


2T{-v) 


cosh(t) : 


- 4 * 


1, — C tanh t. 


+ \[2 7T 


In 2 + ip(—is/ 2) 

2r(-^) 




(56) 


where i^{y) = 4^ lnr(z^) is the Digamma function. 


3. Series expansions of products of parabolic cylinder functions 

We show in this section that, using the integral representation from the previous section, 
we can establish some series expansions of products of parabolic cylinder functions with 
different parameters. 

We use here Tricomi series expansion of confluent hypergeometric function T in terms 
of Laguerre polynomial [16, p. 275, Eq.(3)]. 




1 y 1 {x) 
n + u 


(57) 


Identity (57) is valid for positive real x and for /j, < |. It is also known that Laguerre 
polynomial is a particular case of confluent hypergeometric function T; see [16, p. 268, 
Eq.(36)]. 


(—l) n 

= -j— T(-n,^;x) (58) 

n! 

Appealing now identity (20) for a = 1, representing then \I/(z/, //; x 2 ) in terms 
of Laguerre polynomial as in (57) and finally integrating term by term in the sum, 
we establish the following series expansion of the product of parabolic cylinder functions. 


D_ v (-z)D^ u _i{z) 


1 D n (z)D ll+n -i(z) 
T ^h'o n\(n + y) 


(59) 


where the series of (59) converges for every z£R and for ft < |. 

Transforming (59) with identity (24), it is easy after simplification to find also 


D- v {z)D ^ u ^i(z) 


1 ^ (-lfA t (z)IW-i(z) 

F ( iy )~o n\{n + u) 


(60) 


which is valid for 2 : > 0 and /./ < |. 

If we let n = 1 in (59) and (60), we deduce the series expansion of D_ u (z)D_ u (±z) in 
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terms of Hermite’s polynomials 1 D n (z). 


D_ v (z)D-„(±z) 


1 ) n Dn(z) 2 

F ( iy )~' 0 n\(n + v) 


(61) 


Formula (61) can also be found in Malyshev paper [7] or in its generalization version 
provided by Glasser, in [8], for unrelated arguments y and z. 

Furthermore, from corollaries 1 and 2, it is possible to extract other series expansions of 
the product of parabolic cylinder functions. If we let /j = 1 in (43) and use the generating 
function of Laguerre polynomial L~ l (z) in [5, p. 1002] 


— 7, coth t 
C £ 


2 +°° 

= e-T£L,- 1 (z 2 ) e~ 2nt , 

n =0 


t > 0 


(62) 


we obtain, after integrating term by term and simplifying, the following series expansion 
of D_ u (z) 2 


z 2 +°° n — v—n( 

-%>S^T-ir z 2 W 


D- V (zf = Tre 


n=0 


r ( i ±§±^) 2 


(63) 


Likewise, if we proceed as in (63), we can easily obtain 


, +oo 


D-v{z)D-v(-z) = ire 2 ^(-lfV (z 2 ) 


n=0 


(64) 


Identities (63) and (64) are valid for positive real number z and zz £ C or for all z £ C 
and v is a negative integer. 


4. Conclusions 

In sum, we have proved that Laplace transform of the function x :i-£ 
1 2 

x ,i ~ 1 e 2 X F (y, yL\ x 2 ), where F = $ or T, is explicit and equals the product of parabolic 
cylinder functions with different parameters D u (±z)D u+fJi -i(z). These results lead also to 
Nicholson-type formulas and to new series expansions of the product D v (±z)D v +fj,-i(z). 
In addition, we have established, in theorem 3, Laplace and Fourier transforms of the 
function x :t —y X 2u ^ 1 e~ 2x ~ F (z/, u + \ ;x 2 ), where again F = $ or T. Furthermore, the 
generalization of formulas in theorem 3 seems to do not involve confluent hypergeo¬ 
metric functions. I guess that their generalization would entail higher order confluent 
hypergeometric functions 2 -FT 

Finally, the results of this paper can also be useful for the calculation of integrals involv¬ 
ing the product of two, three or even four parabolic cylinder functions. References 


1 Hermite’s polynomial H n is related to D n by H n (z) = 22 e 2 D n (y/2z). 
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